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ABSTRACT 

We continue the investigation of overcoming Gibbs phenomenon, i.e., obtaining exponen- 
tial accuracy at all points including at the discontinuities themselves, from the knowledge of 
a spectral partial sum of a discontinuous but piecewise analytic function. 

We show that if we are given the first N expansion coefficients of an L 2 function f(x) in 
terms of either the trigonometrical polynomials or the Chebyshev or Legendre polynomials, 
we can construct an exponentially convergent approximation to the point values of 
f(x) in any sub-interval in which it is analytic. 
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1 Introduction 


In this paper we continue our investigation of overcoming Gibbs phenomenon, i.e., recovering 
pointwise exponential accuracy at all points including at the discontinuities themselves, from 
the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function, 
which we started in [5] and [6]. 

Spectral approximations based upon trigonometric polynomials (Fourier, for periodic 
problems) or polynomials (Chebyshev or Legendre, for non-periodic problems) are expo- 
nentially accurate for analytic functions [4], [3]. However, for discontinuous but piecewise 
analytic functions, the spectral partial sum approximates the function poorly throughout 
the domain. Away from the discontinuity only first order accuracy is achieved. Near the 
discontinuity there are 0(1) oscillations which do not decrease with N , the number of terms 
retained in the spectral sum. This is known as Gibbs phenomenon. 

In [5] and [6] we treated a representative problem : 

• Let f(x) be an analytic but nonperiodic function. Suppose that we are given its first 
-N < k < N Fourier coefficients. Can one obtain exponentially convergent pointwise 

approximations in the maximum norm? 

To solve this problem we used the Gegenbauer polynomials C*(x), which are orthogonal 
in [—1, 1] with the weight function (1 — x 2 ) A- ^. 

The procedure of overcoming the Gibbs phenomenon consists of two steps: 

1. Given the Fourier partial sum of the first N terms, we first recover the first 
m ~ N Gegenbauer expansion coefficients with exponential accuracy. This can 
be achieved for any L 2 function, as long as we choose A in the weight function 
of Gegenbauer polynomials to be proportional to N. The error incurred at this 
stage is called the truncation error. 

2. The next step is to prove, for an analytic function , the exponential convergence 
of its Gegenbauer expansion, when the parameter A in the weight function is 
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proportional to the number of terms retained in the expansion. The error at this 
stage is labelled the regularization error. 

Thus we have shown how to overcome the Gibbs phenomenon for an analytic but non- 
periodic function. By a simple shift this procedure covers the case of any analytic function 
with one discontinuity. 

In this paper we treat a broader class of problems. Here we assume that f(x) is an L 2 
function on [—1, 1] and analytic in a subinterval [a, 6] C [—1, 1]. We assume that the spectral 
partial sum (based on either the Fourier or Chebyshev or Legendre expansion) of a function 
over the full interval [—1, 1] is known, and try to recover exponentially accurate point values 
over a subinterval [a, 6]. 

We will follow the same path as in [5]. Basically we will show that the first N Fourier 
(or Chebyshev or Legendre) expansion coefficients contain enough information such that a 
rapidly converging Gegenbauer expansion in the subinterval [a, b] can be constructed. As 
before, we will separate the analysis of the error into two parts: truncation error and reg- 
ularization error. Truncation error measures the difference between the exact Gegenbauer 
coefficients and those obtained by using the spectral partial sum. These will be investigated 
in Section 3 for the Fourier case and in Section 4 for the Chebyshev and Legendre cases. 
Regularization error then measures the difference between the Gegenbauer expansion using 
the first few Gegenbauer coefficients and the function itself in a sub-interval [a, 6], in which 
the function is assumed analytic. The estimation of the regularization error is in Section 5. 
In Section 2 we will provide some useful properties of Gegenbauer polynomials to be used 
later. Section 6 contains a summary theorem and some remarks. 

Throughout this paper, we will use A to denote a generic constant or at most a polynomial 
in the growing parameters. It may not be the same at different locations. 

The message of this paper is that the knowledge of the partial spectra] sum of an L 2 
function in [—1,1] furnishes enough information so that an exponentially convergent approx- 
imation can be constructed in any subinterval in which f(x) is analytic. 
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In particular, if /(x) is piecewise analytic , then an exponentially convergent approxima- 
tion in the maximum norm can be recovered from its partial spectral expansion. 

2 Preliminaries 


This section is devoted to a collection of results about the Gegenbauer polynomials. Though 
the results are classical, they are not widely known and certainly not in our context. We 
rely heavily on the standardization in Bateman [2]. 

We start by defining the Gegenbauer polynomials C*(x) in the following 

Definition 2.1 

The Gegenbauer polynomial C„(x) is defined by 


(i - *y-*ci(») = 



( 2 . 1 ) 


where G( A, n) is given by 


G(A,n) 


+ |)^( ra + 

r(2A)r(n + A + A) 


( 2 . 2 ) 

□ 


Under this definition we have 


and 


C A (1) = 


T(n + 2A) 

n!r(2A) 


|Cj(x)|<C,‘(l), -1<x<1 


(2.3) 


(2.4) 


The Gegenbauer polynomials are orthogonal under their weight function, in fact we have 

Lemma 2.2 

The Gegenbauer polynomials satisfy the following orthogonality condition: 

fj 1 - x ! ) A -‘ c?(x)c;(z)ix = h.„h x n (2.5) 


where 




r(A + A) 

r(A)(n + A) 


(2.6) 
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□ 

We will need to use heavily the asymptotics of the Gegenbauer polynomials for large n 
and A. For this we need 
Lemma 2.3 (Stirling) 

For any number x such that x > 1 we have 

(2n)^x x+ ^e~ x < r(a; + 1) < (2 tt)» x x+ i e~ x (2.7) 

□ 


Lemma 2.4 

There exists a constant A independent of A and n such that 


i-i 


c. A (i) < < AjAL^c^) 


(n + A) 


(n + A) 


The proof follows from (2.6) and Stirling’s formula (2.7). 


( 2 . 8 ) 


□ 


In the analysis we will extensively use the relationship between the Fourier functions 
e iknx and the Gegenbauer polynomials C*(x). In particular, we will focus our attention on 


the Legendre polynomials 




P n {x) = ck{x) 

(2.9) 

and the Chebyshev polynomials 




Tn{x) = \C° n {x) 

(2.10) 


Our basic formula is taken from [2], page 213: 


Lemma 2.5 


1 

hf 


J\l - x 2 ) x ~^e invx Cl'{x)dx = T(A) i‘(l + A)J, +A (7rn) 


where J^x) is the Bessel function. 


( 2 . 11 ) 
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□ 


Lemma 2.5 gives us the Gegenbauer expansion coefficients of the Fourier functions. How- 
ever, what we really want are the Fourier coefficients of the Gegenbauer polynomials. We 
have those coefficients, luckily, for the Legendre and Chebyshev polynomials. 

Lemma 2.6 

Let a jf be the Fourier coefficients of the Legendre polynomial Pn(x ) defined in (2.9) and 
be the Fourier coefficients of the Chebyshev polynomial Tpj(x) defined in (2.10), i.e. 


and 


Then 


and 


OO 


«*(*)= £ <*?«“" 
k— — oo 

(2.12) 

OO 

r»(*)= £ 

k= — oo 

(2.13) 

i N 


= VTk JN ^ k) 

(2.14) 

h = 2 ^ JV ( jr *) 

(2.15) 


□ 


Equation (2.14) follows from (2.11) upon substituting A = | and using (2.6) and (2.3). 
Equation (2.15) can be found in [7], page 836. 

□ 

Lemma 2.6 provides us with a better understanding of the terms “small scales” and 
“large scales”. The question is the following: it is clear that e ,w ” is small scale for large JV, 
intuitively T^(x) and Pn(x) are also small scale functions, however the expansions (2.12) 
and (2.13) contain all the scales. So, in what sense are Tpj(x) and Pn(x) small scales? 

The answer is given in the following Lemma: 

Lemma 2.7 
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Let and be the Fourier coefficients of the Legendre polynomial P/v(x) and the 
Chebyshev polynomial Tjv(x) respectively, given in (2.12)-(2.15). Then 


• (t ( evk \^ 


|aj^ | < ylmin 
|fef| < A min 

where A is a constant independent of k or N. 

Proof: 


(2.16) 


(2.17) 


We start with the following estimate for the Bessel function J n (nz ) (see [1], page 362): 

(2.18) 


|j„(nz)| £ 


ze 


\/T — 2 


l +VT= 


0< z <1 


From (2,18) we can deduce trivially that 


f ez x 


P.MI < (y)“ 0 < z < 1 


(2.19) 


We also have the uniform bound on Bessel functions 


\Jn(z)\ < 1 (2.20) 

Using (2.19) and (2.20) in (2.14) and (2.15), with z = ^ and n = N, we obtain (2.16) 
and (2.17). 

□ 

Lemma 2.7 provides the answer why Pn(x) and Tn(x) (for large N) can be considered 
as small scale functions. Basically only the terms that are greater than N appear in their 
Fourier expansions (2.12) and (2.13). The lower terms decay exponentially with N. 
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3 Truncation Error in a Sub-interval — Fourier Method 


Consider an arbitrary L 2 function f(x) defined in [—1,1]. Suppose that the first 2N + 1 
Fourier coefficients of this function are given 

/(*) = \ J\ f(x)c~ ik ™dx |A:i<7V (3.1) 

In [5] we showed that we can recover (with exponential accuracy) the first m ~ N terms 
in the Gegenbauer expansion of f(x) based on the interval [—1,1]. 

Here we are interested to find the Gegenbauer expansion of f(x) based on a sub-interval 
[a,b] C [-1,1]. This will allow us to handle multiple discontinuities and discontinuities of 
unknown locations. We will show that in this case too, we can get exponential accuracy. We 
start by introducing the local variable 
Definition 3.1 

The local variable £ is defined by 


x = z(0 = + S 


(3.2) 


where 



Thus when a < x < b, — 1 < £ < 1. 



(3.3) 


□ 


Denote now the Fourier partial sum by //v(x), namely 

/»(*)= E /(t)e'*’ x (3.4) 

k=-N 

As mentioned before, we are trying to recover the Gegenbauer expansion coefficients 
based on the sub-interval [a, b], i.e. we would like to find the first m ~ N coefficients /*(/) 
in the expansion 

CO 

/MO) = £/,*(<><?, >(0 (3.5) 

/=o 
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where 


/, A (0 = ^ jj 1 - ( a ) A -*/(x(0)C?«)de (3.6) 

Of course we do not have /<*(/) at our disposal but only an approximation based on the 
Fourier partial sum //v(x), thus we have 

m = A /_',(! -fV'bwWOJCftfK (3.7) 

The Truncation Error is the measure of deviation of the approximate Gegenbauer expan- 
sion from the true expansion: 

m 

TE(\,m,N,e)= max | £(/*(0 - ft A (OK?(OI (3.8) 

(=0 

We can show that the truncation error is exponentially small: 

Theorem 3.2: 

Let f(x ) be an L? function whose first 2N 1 Fourier coefficients f(k) defined in (3.1) are 
known. Let g^(l) be the Gegenbauer coefficients, based on the sub-interval [a, 6] expansion 
of the Fourier partial sum /)v(x) defined in (3.7), and f c (l) be those of f(x) defined in (3.6). 
Assume that 

A = aeN m = fieN (3-9) 


Where e is the half length of the sub-interval [a, 6] defined in (3.3). Then the truncation 
error is bounded by 

TE(atN, /3tN, N, e) < Aq cN , (3.10) 


where 


9 = 


(P + 2a)^ 


2a 


( 2 , Ke) a a a fiP 

and A grows at most quadratically with N . In particular, if a = ft « L then 


(3.11) 


q ~ 0.8 


Proof: 



The proof follows exactly the one in [5], Theorem 3.1 and 3.3, with proper accounting 
for the interval length e. 

□ 

4 Truncation Error in a Sub- interval — Chebyshev and 
Legendre Methods 

In this section we assume that we have the first N + 1 Chebyshev or Legendre coefficients 
in the expansion of an L 2 function f(x) based on the interval — 1 < x < 1. We show how 
to get its Gegenbauer expansion based on a sub-interval [a, 6] C [—1,1] with exponential 
accuracy in N. The proof makes use of Lemma 2.7, which states that the low modes of the 
Fourier expansion of the Chebyshev or Legendre polynomial of degree N decay exponentially 
with N. Thus only the high modes count, and those are shown in Section 2 to have small 
Gegenbauer coefficients. 

For simplicity of notation we will deal with C k (x), which is the Legendre polynomial 
P k (x) (see (2.9)), and C k (x), which is - times the Chebyshev polynomial T k (x) (see (2.10)). 
The Chebyshev or Legendre expansion coefficients of f(x) are given by 

/"(*) = T j\ 1 - x 2 r-lf(x)C£(x)dx (4.1) 

for n = 0 or fi = 

In this section we consider functions f(x) satisfying 

Assumption 4.1: 

|/ M (A;)| — A for all k. 

□ 

We remark that if f(x) is an L 2 function the assumption is fulfilled. 

We assume that we know the first N + 1 Chebyshev or Legendre coefficients, /**(&) for 
0 < k < N, and define 

«(*)=£/'(*)<?(*> ( 4 - 2 > 

k=0 
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Note that f^/(x) does not converge fast to f(x ) if there exist discontinuities inside the domain. 

The function f(x) has also a Gegenbauer expansion in a sub-interval [a, 6]. With £, e and 
8 defined in (3.2)-(3.3), we have 

OO 

/w + «) = £/ A (')c A (e) (4.3) 

1=0 

where the Gegenbauer coefficients /*(/) are defined by 

;, A (o = A J\ i - (4.4) 

As before, we do not have / A (/) at our disposal, but only an approximation based on the 
partial Chebyshev or Legendre sum /^(x), thus we have 

s, A (o = ujj 1 - + <)c, A (Oo(e (4.5) 

How well do <7 A (/) approximate / f A (/)? To answer this question we define 

Definition 4.2: 

The truncation error is defined by 

m 

T E(X, m, N, e) — max | ^(/ A (/) - ^(0)^(0! (4-6) 

- cs 1=0 

where / A (/) are defined by (4.4) and <) A (/) are defined by (4.5). 

□ 

The truncation error is the measure of the distance between the true Gegenbauer expan- 
sion in the interval [a, 6] and its approximation based on the Chebyshev or Legendre partial 
sum in [ — 1, 1]. 

We first have the following lemma: 

Lemma 4.3: 


TE(X,m, N,e) < £ \f>{q )\ £ 

q = N + 1 1=0 




(4.7) 
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Proof: 


From (4.4) and (4.5) we have 


|/. A (') - », A (o| = /ji - e?-Hm + *) - + ml 


Substituting (4.8) into (4.6), recalling (2.4) and 


m+s) -/&{'(+*)= £ rwc'w + f) 

q=N + \ 


we obtain (4.7). 


For simplicity of notations we denote: 

|C, A (1) 


Cq.fA = 


q> N 


(4.8) 


(4.9) 


□ 


(4.10) 


hT , 

In order to estimate this term, we consider the Fourier series approximation of C^*(e£ + 6) 
from (2.12) and (2.13) to get 

Lemma 4.4: 


OO / O \ A 

F , < /lC, A (l)r(A) p )— ) |J,+a(*x£)| Kl (4.11) 

where a k is given either by (2.14) or by (2.15), and A grows at most as fast as a linear 
polynomial in A and /. 

Proof: 

We expand Cjf(x) using (2.12) or (2.13): 


C?(x) = jr a q ^e ,krx = f; a£V fc "V tarf 

k— — oo k ~ — OO 

where a k 2 = a q k is given by (2.14) and a q k ° = ^b q k is given by (2.15). 
We now use (2.11) to get 

1 ,i 1 / ^ xA 


(4.12) 


E/hl -{ I ) 1 -ie i ‘"<e M C l , ((W =r(A)(]j^) (1 + A) |J, + »(fcr«)| (4.13) 


Substituting (4.12) and (4.13) into (4.10) gives us (4.11). 


11 



□ 


We are ready now to estimate directly F q ^\ defined in (4.10). 

Lemma 4.5; 

Consider Fgj^ defined in (4.10). Let A = and l = jX where 7 < 1 . Then for every 


Fni.x < Ae x 




(4.14) 


q > N we have 

\,/,a ^ 

where A grows at most quadratically with N. 

Proof: 

Let k q = be the integer part of We split the sum on the right hand side of (4.11) 
into 

2 ' A 


^,a<^(E+ e 

k=k<i+l) 


c,A(1)r(A) (lL) Wl = A (f£s + F $*) ( 4A5 ) 


We first use (2.3) and (2.7) to obtain: 


7 ( 1 ) < A 


'(7 + 2 . 727 


47 "! 




and 


T(A) < AA A e 


A„— A 


(4.16) 


(4.17) 


We start with an estimation of the second sum 

CT(1)1(A)( 1 2 v 

k=k q +l 

We use (4.16), (4.17), the fact that |«/j + >,(:r)| < 1 and |aj[| < 1, and k > k q = 


f$„ = E 7(l)r(A)(^) A |J,«(^e)||a 


and 


A = to get 




(4.18) 


where A grows at most linearly with N, the growth coming from the summation. For the 
first sum 


7% = E 7(W) (j^) A w«(A«)| 1 a 


21. 
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we still use (4.16), (4.17) and the bound \J t+ \(x)\ < 1. However, for the estimate on \a q k \, 
we use (2.16), (2.17): 


Kl < , (f )' 


and the fact q > N > A to get 


and, upon substituting for k q = 


i (enk\< ( 27A V ( e*k,V ( 27A V 

2enkeJ q \ 2q ) \2e7r k q e ) 


and A = 






Combining (4.18) and (4.21) we obtain (4.14). 


(4.19) 


(4.20) 


(4.21) 


□ 


We are now ready for the main theorem of this section: 

Theorem 4.6: 

Let the truncation error be defined in (4.6). Let A = and m < A, then 

TE(\,m,N,t) < Ae~ x (4.22) 

where A grows at most as a fourth order polynomial in N. 

Proof: 

The theorem follows from (4.7), the Assumption 4.1, and (4.14). 

□ 


5 Regularization Error for a Sub-interval 

In this section we study the second part of the error, which is caused by using a finite 
Gegenbauer expansion based on a sub-interval [a, 6] C [—1,1], to approximate a function 
f(x) which is assumed analytic in this sub-interval [a,f>]. Since in previous sections we have 
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established exponentially small truncation errors if A and m are both growing linearly with 
N, in this section we will consider the case of A ~ m. 

We will assume that f(x) is an analytic function on [a, 6] satisfying 

Assumption 5.1 

There exists constants p > 1 and C(p) such that, for every k > 0, 

<C(p)- k (5.1) 
P K 

□ 


max 

a<r<b 


d k f 


dx k 


(x) 


This is a standard assumption for analytic functions, p is the distance from [a, b] to the 
nearest singularity of f(x) in the complex plane (see, for example, [8]). This assumption and 
the following proof can be modified to allow for an arbitrary ellipse region around [a, b] for 
the analyticity of f(x) (H. Vandeven, private communication). 

Let us consider the Gegenbauer partial sum of the first m terms for the function f(e£ + 8): 


/*'■«) = ££(')<?«) 


(5.2) 


1=0 


with e and 6 defined by (3.2) and (3.3), and the Gegenbauer coefficients based on [a,b] 
defined by 

/*w= p /_'(i - e y-*m+ <)c, A «Me (5.3) 


We want to estimate the regularization error in the maximum norm: 


RE(\,m,e)= max 


/« + «) - E/. a (')c, a K) 

/=0 


(5.4) 


We have the following result for the estimation of the regularization error, when A ~ m: 


Theorem 5.2; 

Assume A = 7 m where 7 is a positive constant. If f(x ) is analytic in [a, 6] C [—1, 1] 
satisfying the Assumption 5.1, then the regularization error defined in (5.4) can be bounded 

by 

RE(-fm,m,e) < Aq m (5.5) 
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where q is given by 


which is always less than 1 . 
constant, then 


9 = 


e(l + 27 ) 1+2 ' y 


(5.6) 


pi'* 

In particular, if 7 = 1 and m = /37V where /? is a positive 


RE({3N, f3N, e) < A?* 


(5.7) 


with 



(5.8) 


Proof: 

The proof follows in exactly in the same fashion as in [5], Lemma 4.2, Theorems 4.3 and 
4.4, with proper accounting for the transformation x(£) = + 8. 


□ 


6 Concluding Remarks 


In this Section, we combine the estimates for truncation errors and regularization errors in 
previous sections, to obtain the main theorem of this paper: 

Theorem 6.1: Removal of the Gibbs Phenomenon for the sub-interval case 

Consider an L 2 function f(x) on [—1, 1], which is analytic in a sub-interval [a, b] C [—1, 1] 
and satisfies 


max 

a<x<b 


d k f 


dx k 


(*) 




Assume that either the first 2N -f 1 Fourier coefficients 

/(*) = \ j\ f(x)e~ ik **dx , -N < k < N 
or the first N + 1 Chebyshev or Legendre coefficients 

nk) = ^ fj 1 - x'r-lfWCZWdx, ti = 0 or IX = J 


( 6 . 1 ) 
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are known. Let 0 < / < m be the Gegenbauer expansion coefficients, based on the sub- 

interval [a, 6], of the spectral partial sum /n(x) in (3.4) or /£(x) in (4.2). These Gegenbauer 
coefficients are defined in (3.7) for the Fourier case and in (4.5) for the Chebyshev or Legendre 
case. Then for A = m = fitN where /? < ^ (Fourier case) or /? < ^ (Chebyshev or 
Legendre case), we have 


max 
— 1 <£<1 


m 

M + Q-EtfM 1 K) 


/=0 


< A + q^j 


where 



for the Fourier case and 



( 6 . 2 ) 


for the Chebyshev or Legendre case, and A is at most a fourth order polynomial in N. 


Proof: 

Just combine the results of Theorems 3.2, 4.6 and 5.2. 


□ 


We finally give two remarks: 

Remark 6.2 

It can be seen from (6.2) that the truncation error gets bigger when the interval half- 
length e decreases, because the parameters m and A must decrease with t. This in turn also 
affects the regularization error, although there is a factor e in qn which offsets this effect. 

□ 


Remark 6.3 

No attempt has been made in this paper to optimize the parameters. 


□ 
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